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We have investigated the thermal equation of state of bcc tantalum from first principles us- 
ing the full-potential Linearized Augmented Plane Wave (LAPW) and mixed-basis pseudopotential 
methods for pressures up to 300 GPa and temperatures up to 10000 K. The equation of state at 
zero temperature was computed using LAPW. For finite temperatures, mixed basis pseudopotential 
computations were performed for 54 atom supercells. The vibrational contributions were obtained 
by computing the partition function using the particle in a cell model, and the the finite temperature 
electronic free energy was obtained from the LAPW band structures. We discuss the behavior of 
thermal equation of state parameters such as the Griineisen parameter 7, the thermal expansivity 
a, the Anderson-Griineisen parameter St as functions of pressure and temperature. The calcu- 
lated Hugoniot shows excellent agreement with shock-wave experiments. An electronic topological 
transition was found at approximately 200 GPa. 



I. INTRODUCTION 



We investigate from first-principles the thermal equa- 
tion of state of body-centered cubic (bcc) tantalum, a 
group V transition metal, which is a useful high-pressure 
standard due to its high structural mechanical, thermal 
and chemical stability. Ta has a very high melting tem- 
perature, 3269 K at ambient pressure, and its bcc struc- 
ture is stable for a large pressure range. Static diamond 
anvil cell experimental up to 174 GPa and full-potential 
linearized muffin tin orbital (LMTO) calculations! up to 
1 TPa conclude bcc phase of Ta is stable at these pres=. 
sure range. Similarly, shock compression experiments^! 
showed melting at around 300 GPa, but no solid-solid 
phase transition. 



II. STATIC EQUATION OF STATE 

Firstly, we discuss the static high pressure properties 
of Ta, which we obtained from first principles by using 
the LAPW methodoB. The 5p,4f,5d and 6s states were 
treated as band states, others are described as core elec- 
trons. n We used both the local density approximation 
(LDA)tJ and generalized gradient approximation (GGA)EI 
for the exchange-correlation .potential. The Monkhorst- 
Pack special fc-point schemea with a 16x16x16 k-point 
mesh (140 k-points within the irreducible Brillouin zone 
of the bcc lattice) was used after convergence tests. The 
convergence parameter KK max was 9.0, and the muffin 
tin radii were 2.0 bohr, giving about 1800 plane-waves 
and 200 basis functions per atom at zero pressure. The 
total energy was computed for 20 different volumes from 
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62.5 to 164 bohr 3 (1 bohr=0. 5l29177 A), and the energies 
were fit to the Vinet equation^, 

E(V,T) = E (T) + 9R o( T )MT) {1 + _ x) _ 1} (1) 

xexp{t; (1 - x)}} 

where Eq and Vq are the zero pressure equilibrium energy 
and volume respectively, x = (^r-) 3 and £ = ^(K' — 1) , 
K (T) is the bulk modulus and °K' Q {T) = {dK(T)/dP) . 
The subscript alone throughout represents the stan- 
dard state P=0. All equations of state here are for an 
isotherm or static (T=0) conditions, unless specified oth- 
erwise. Pressures were obtained analytically from 

P(T, V) = { 3^o(D(l-x) _ 
x z 

The calculated equation of state is compared with exper- 
iments!^ in Fig. 0. The LAPW GGA results are found 
to be more accurate than the LDA. The discrepancies are 
larger beween theory and experiment at high pressures; 
this may be due to strength effects in the experimentsEil, 
since we find good agreement with the experimental 
Hugoniot to 400 GPa (discussed belowl-. For the pseu- 
dopotential mixed-basis calculationsE3ll3 (MBPP, dis- 
cussed below), we find that the LDA agrees fairly well 
with the experiments, indicating compensating errors be- 
tween the pseudopotential and the LDA. Since we are 
using a first-principles approach, and want to avoid ad 
hoc variations in procedure to get better agreement with 
experiments, we use the computational more accurate 
method with the least approximations, that is LAPW 
and GGA, rather than the MBPP with LDA, in spite of 
the fortuitously better agreement of the latter with the 
room temperature data. 

The residuals between the calculated and fitted ener- 
gies show large deviations for volumes less than 80 bohr 3 
(Fig. ||). Note that the residuals are all small and not ev- 
ident for the large energy scale shown in Fig. (?]a. Other 
equation of state formulations, such— as the extended 
Birch equationli-3 show the same trendt^l. When the fit is 
restricted to volumes greater than 80 bohr 3 , or extend- 
ing the Vinet equation by two more parameters related 
to the next two pressure derivative of bulk modulus, it is 
improved significantly. Hence, these large residuals are 
related to different high- and low-pressure behavior of Ta. 
The band structures and densities of states (Fig. ||) show 
a major reconfiguration of the Fermi surface. This elec- 
tronic topological transition is the reason for this change 
in compression, and the behavior of the residuals for the 
fitted equations of state. This indicates that systematic 
deviations from simple equations of state can be used to 
find subtle phase transitional^. 

Spin-orbit interactions may be important for Ta. In 
order to test this, we included spin-orbit coupling by sec- 
ond variational treatments including 20-80 -bands. In 
contrast to the fully relativistic LMTO resultsB, we found 
only negligible effect on the equation of state, so our com- 
putations were done without spin-orbit. 



III. CALCULATION OF THERMAL PROPERTIES 



In order to compute the high temperature properties 
of Ta, we separated the Hclmholtz free energy aaL3: 

F{V,T) = E static (V) + F el (V,T) + F vib (V,T), (3) 

where E sta tic{V) is the static zero temperature energy, 
F e i(V,T) is the thermal free energy from electronic exci- 
tations, and Fvib{V,T) is the vibrational contribution to 
the free energy. E sta tic{V) and F e i(V, T) were computed 
using the LAPW method with the GGA. F vib (V : T) was 
computed using the particle-in-a-cell (PIC) model with a 
mixed basis pseudopotential method, as described below. 
The main differences we find between LDA and GGA are 
in the energy versus compression, but we find only small 
differences in energy versus atomic displacements. Since 
LDA converges much faster than GGA, and is also faster 
per iteration cycle, we used LDA for the large supercell 
computations required for the vibrational contributions. 
Differences from using GGA for this are negligible. 
The electronic part F e i(V, T) of the free energy is: 

F el {V, T) = E el {V, T) - TS el {V, T) (4) 

where E e i{V,T) is the internal energy due to thermal 
electronic excitations, 

S a (V, T) = -2k B fi ^ fi + (1 - fi) In(l - /*) (5) 



is the electronic entropy, and the Fermi-Dirac occupation 
fi is 

ti are the eigenvalues, fi is the chemical potential, and ks 
is the Boltzmann constant. The vibrational free energy 
is given in terms of the partition function as 

F vib (V,T) = -k B TlnZ. (7) 



The particle-in-a-cell modet-irEJ was used to calculate 
the partition function. In this model, the partition func- 
tion is factored by neglecting atomic correlations. An 
atom is displaced in its Wigner-Seitz cell in the poten- 
tial field of all the other atoms fixed at their equilibrium 
positions, i.e. the ideal, static lattice except for the wan- 
derer atom. The partition function is simply a product 
of identical functions for all the atoms, involving an in- 
tegral of Boltzmann factor over the position of a single 
atom inside the Wigner-Seitz cell, 



Z C eii = \~ 3N { I exp 



(U(r) - Uo(T )) 



dfr, (8) 



where A = h/i^nmkBT) 1 ^ 2 is the de Broglie wavelengths 
of atoms and U(r) is the potential energy of the system 
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with the wanderer atom displaced by radius vector r from 
its equilibrium position. The advantage of the cell model 
over lattice dynamics based on the quasiharmonic ap- 
proximation is that anharmonic contributions from the 
potential-energy of the system have been included ex- 
actly without a perturbation expansion. On the other 
hand, since the interatomic correlations between the mo- 
tions of different atoms is ignored, it is only valid at tem- 
peratures above the Debye temperature. Diffusion and 
vacancy formation are also ignored, so premelting effects 
are not included. We have used the classical partition 
function, so quantum phonon effects are not included. 
Thus the heat capacity and thermal expansivity do not 
vanish at low temperatures, for example. The present 
results are appropriate for temperatures above the De- 
bye temperature (245 K in Tata) and below premelting 
effects. Since the Debye temperature is below room tem- 
perature in Ta, the classical thermal properties should be 
reasonable even down to room temperature. 

The electronic thermal free energy was obtained using 
the Mermin theorems (Eq. ||). The charge density is 
temperature dependent through both occupation num- 
bers according to the Fermi-Dirac distribution and self- 
consistency. The electronic contributions to the thermal 
free energy were computed by the LAPW method us- 
ing the same computational parameters as the T = OK 
computations described in Section |n]. 

For the vibrational contributions, it is necessary to do 
a large number of large supercell calculations, which is 
computationally intractable by the LAPW method, but 
is achievable with the MBPP methodic. In this mixed- 
basis approach, pseudo-atomic orbitals and a few low- 
energy plane waves are used as the basis set within a 
density functional, pseudopotential calculation. It was 
shown that the method offers a computationally efficient 
but accurate alternative. 

A semi-relativis.tic, nonlocal and norm-conserving 
Troullier-MartinsEil pseudopotential (with associated 
pseudo-atomic orbitals) was used to describe the Ta 
atoms. The pseudopotential was generated from an 
5d 3 6s 2 6p° atomic configuration with cutoff radii 1.46, 2.6 
and 3.4 bohr for 5d, 6s, and 6p potentials, respectively, 
with non- linear core corrections. The cutoff radii were 
optimized by testing the transferability of the pseudopo- 
tential by considering the reasonable variations of the 
reference atomic configuration and by comparing the log- 
arithmic derivative of the pseudo-wavefunctions with all- 
electron values in the valence energy range. The 6s po- 
tential was chosen as the local component while 5d and 6p 
were kept as nonlocal while transforming the potential to 
the nonlocal separable Kleinman-Bylander formca. A full 
plane-wave representation is used for the charge density 
and potential, and a smaller cutoff is used for the basis 
set. The pseudoatomic orbitals are expanded in the large 
plane wave set for evaluation of the potential and charge 
density integrals in the Hamiltonian and overlap matri- 
ces, and in the total energies. After checking the energy 
convergence, 550 eV and 60 eV were used for the large 
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and small energy cutoffs, respectively, in the solid cal- 
culations. The exchange-correlation effects of electrons 
were treated within LDA. The T = OK equation of state 
of bcc Ta was computed to test the pseudopotential, and 
is compared with experiment and LAPW results in Fig. 

For the PIC computations, a supercell with 54 atoms 
was used. The MBPP calculations were carried out 
on this 54 atoms supercell using LDA for exchange- 
correlations effects and 4 special k points for BZ inte- 
grations. The potential energy surface was then calcu- 
lated as a function of the displacements of the wanderer 
atom. Symmetry was taken into account in order to re- 
duce the number of computations. The integrand in Eq. ^ 
has a Gaussian-like shape and decays rapidly, and essen- 
tially is zero at half of the interatomic distances even 
at very high temperatures. Therefore, integration over 
the Wigner-Seitz cell can be replaced by an integration 
over the inscribed sphere. Also, the radial part of the 
integrand is invariant under point group operations of 
the lattice, hence a numerical quadrature can be used 
for angul ar .. in tegration based on the method of special 
directionsOO. In this method, the radial integral is ex- 
panded in terms of lattice harmonics, cubic harmonics 
for a cubic lattice, then a quadrature rule is derived for 
the angular integration in terms of the radial integration 
by choosing special directions, r^, in such a way that the 
contribution from I =/= terms, as many lattice harmonics 
as possible, is zero. In all of the computations, we used 
one special direction which integrates exactly up to I = 6 
cubic harmonicsEj. Then, the potential energy was cal- 
culated at 4-6 different displacements along this special 
direction. In order to model the potential, these com- 
puted values were fit to an even polynomial up to order 
8, which shows the anharmonicity very clearly, since a 
second order fit describes the data poorly. Finally, the 
cell-model partition was calculated from Eq. |§| by carry- 
ing out the integration numerically, and the vibrational 
free energy is simply given in terms of partition function 
by Eq. 0. 

IV. THERMAL EQUATION OF STATE 

We treated the resulting free energies^ three differ- 
ent ways. Firstly, F-V isotherms were fit using the Vinet 
equation of state, giving E (T), V (T), K Q (J) and K' {T) 
as the parameters of the fit (Table 0, Fig. Q). Because of 
the thermal expansivity, the minimum energy shifts to 
higher volumes with increasing temperature, and above 
T=6000 K the P=0 volume is not in the range of volumes 
we studied. The experimental P=0 melting temperature 
is 3270 K, so temperatures above this are non-physical 
in any case. The parameters for higher temperatures are 
Active parameters that describe the higher pressure equa- 
tion of state accurately. The Vinet parameters could be 
fit to polynomials versus temperature to obtain a thermal 
equation of state, but the following approaches require 
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fewer parameters. 

A second thermal equation of state was obtained by 
analyzing the thermal pressure obtained from the Vinet 
fits (i.e. the differences in pressures between isotherms). 
The thermal pressure as a function of volume and tem- 
perature is shown in Fig. ||. The volume dependence 
of the thermal pressure for Ta is very weak up to 80% 
compression. Deviations at higher compressions may be 
due to the fit through the electronic topological transi- 
tion discussed above, and thus due to inflexibility in the 
Vinet equation, rather than a real rise in thermal pres- 
sure. The thermal pressure is also quite linear in T. The 
thermal pressure changes are given by 

P(V,T)-P(V,T )= [ aK T dT (9) 
J To 

and olKt is quite constant for many materials^ in the 
classical regime (above the Debye temperature). Hence, 
a greatly simplified thermal equation of state is to add a 
thermal pressure 

P th {T) = aT (10) 

to the static pressure P static- Since our equation of 
state is classical, we can use this expression at all tem- 
peratures. This simple equation of state gives q = 

( oinv ) = 1 anc ^ (t5t")v = 0' which therefore are good 
approximations over this pressure and temperature range 
for Ta. 

The thermal pressure was averaged over volumes from 
60 bohr 3 to 220 bohr 3 , and is shown as a function of 
temperature in Fig. ||. The solid line has a slope a 
of 0.00442 GPa/K. This is close to aK 0T (T), which is 
0.00460 GPa/K at 1000 K and zero pressure. So a sim- 
ple equation of state for Ta is the static pressure given 
by V (T=0)=123.632 bohr 3 , K (T=0)=190.95 GPa, and 
Kf,(T=0)=3.98 in the Vinet equation (Eq. |) plus the 
thermal pressure P t h = 0.00442T. 

Thirdly, an accurate high temperature global equation 
of state was formed from the T = OK Vinet isotherm and 
a volume dependent thermal free energy F t h as: 

«=3,j'=3 

F th = J2 A ij Ti V j -3k B TkiT (11) 

i=l,j=0 

This is the thermal Helmholtz free energy per atom, 
which must be added to equation ^ to obtain the to- 
tal free energy. The parameters Aij are given in SI units 
in Table ||, which gives the free energy in Joules/atom; 
T is in Kelvin and V in m 3 /atom. The term in TlnT 
is necessary to give the proper classical behavior at low 
temperatures, since we are evaluating the classical parti- 
tion function, with Cy = 3/cb and S = -co at T = OK. 
For the best overall accuracy, the T=0 isotherm was also 
included in the global fit. The global fit is compared with 
the computed free energies in Fig. [?]. The r.m.s. devia- 
tion of the fit is 0.4 mRy. At low temperatures (0 and 
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1000 K) the residuals (Fig. |b) are larger due to the elec- 
tronic topological transition discussed above (Fig. |J) , but 
at higher temperatures this anomaly is less pronounced 
due to the thermal smearing, and the equation of state 
fits well. 

Thermal equation of state parameters such as 
P,a,j,dT,Q, and the heat capacity Cy and Cp can be 
obtained from the global fit by differentiation and alge- 
braic manipulation (seeEil for a collection of useful formu- 
las.) We now discuss the behavior of these parameters. 
The thermal expansion coefficient is presented, and com- 
pared with zero pressure experiments^ in Fig. ||. The 
deviations at lower temperatures are due to the use of 
the classical partition function. The thermal expansivity 
is a quite sensitive parameter, and the errors at mod- 
erate temperatures, which are typical in first-principles 
computations, may come from a number of sources (er- 
ror in the P=0 volume, LDA, the pseudopotential, the 
PIC model, or convergence in k-points or basis set). The 
divergence in behavior at higher temperatures is not due 
to vacancy formation, since the vacancy formation in Ta 
is high (3.2 eV)E3, and the fraction of vacancies at the 
melting point is less than 10 -4 . The temperature range 
over which the anomaly occurs seems too large to be a 
premelting effect. One possible explanation would be an 
incipient solid-solid phase transition in Ta, which would 
not be detected in the PIC method. The upturn in a 
with increasing temperature is apparently a low pressure 
phenomenon. It is possible that it is due to an experi- 
mental problem, such as oxidation of the sample. 

The thermal expansivity drops rapidly with increas- 
ing pressure (Fig. |§|), and this is parametrized by the 
Anderson- Griineisen parameter (Fig. ^ 

d In a 



(12) 



The behavior of St is complex. At low pressures it 
increases with increasing temperature, but at elevated 
pressures it decreases withjtemperature. The parame- 
ter St can be fit to a formEj St = <$t(?7 = 1)?? K where 
V = V/V {T ). The average S T (averaged from 0-6000 K) 
decreases with compression, and a power law fit gives 
S T {jl) = 4.567/ 1 - 29 ; at 1000 K St — 4.757y 117 . Msg-, 
estingly these values are not that different from MgOea 
(S T (v = 1, 1000-ff)) = 5.00 and k = 1.48. The behavior 
is much different than for Fe, where St is constant to 
150 GPa with values of 5.2 and 5.0 for fee and hep, re- 
spectively, after which it drops more slowly than a power 
lawtZl. The difference between St and K' is an important 
anharmonic parameter, and is related to the change in 
the bulk modulus with temperature at constant volume 
and the thermal pressure with compression at constant 
temperature: 

S T -K'=( 91 :^) (13) 



(14) 



o,Kt 



SlnV 
1 fdK T \ 
\~dTJ 



7 



Fig. ||d shows that St — K' is quite small over a large tem- 
perature and compression range, but increases at high T 
and P. This is consistent with the accuracy of the sim- 
ple equation of state (Eq. |Io| ). The behavior of St — K' 
is also surprisingly similar to the behavior of MgO (see 
Fig. 3.3 in||). 

Changes in thermal pressure Pth are given by aKr = 
(§?)y which is shown in Fig. |ic]. Changes in aKr are 
small, but it is interesting that the sign of the change 
with temperature is strongly dependent on pressure, in- 
dicating that experiments to determine ( d T ) at low 
pressures may not be applicable to a very large pressure 
range. 

A most important parameter, particularly for reduc- 
tion of shock data, is the Gruneisen parameter 

v (dP\ aK T V 

^ v \m) v = — (15) 

where E is the internal energy. The Gruneisen parame- 
ter is used in the Mie- Gruneisen equation of state, which 
assumes 7 independent of temperature. Then the ther- 
mal pressure on the Hugoniot, for example is given by 
the change in internal energy by 

Phug Pstatic — "^"(-^Tmg -^siaiic)- (-^) 

Fig. [ll] shows that at elevated pressures, 7 is moderately 
temperature dependent, and it varies more strongly with 
temperature below 100 GPa. The variation of 7 with 
pressure is given by 

dhi7 ri »s 
q = MnV 

which is shown in Fig. [l2| The parameter q is not con- 
stant, as is often assumed, but decreases significantly 
with pressure and temperature. If aKx and Cy were 
constant, Eq. [l5] shows that 5 = 1. Fig. [To] shows that 
aKx is quite constant, so that large changes in q must 
be due primarily to changes in the heat capacity Cy. 

Fig. [13| shows indeed that the heat capacity is a strong 
function of temperature and pressure. This is due mainly 
to the electronic contributions. The experimental Cp at 
zero pressured is also shown. Other than the large differ- 
ences from experiment at very low temperatures, due to 
neglect of quantum phonon effects in the present model, 
there is a large increase in the experimental heat capac- 
ity with increasing temperature that is not seen in the 
PIC results. A similar large increase in the experimen- 
tal thermal expansivity is not predicted by the model 
(Fig. ||). Vacancy formation (not included in the PIC 
model) seems an unlikely source for this discrepancy as 
discussed above. An incipient phase transition or sample 
oxidation seems the most likely cause of the observed be- 
havior in the thermal expansivity and the heat capacity. 

To compare with experiment at high pressures 
and temperatures, we consider the high temperature, 
high pressure equation of state obtained by shock 



8 



compression^. The pressures, Ph and temperatures, Tjj, 
on the Hugoniot of Ta are given by the Rankine-Hugoniot 
equation: 

^Ph(V o (T o )-V) = Eh-E o (T = 0). (18) 

We solved the Rankine-Hugoniot equation using our 
equation of state results by varying the temperature at a 
given volume until it is satisfied. The calculated Hugo- 
niot shows very good agreement with experimental data 
as seen in Fig. |l4]. 

Computations of the Hugoniot using a modified free.- 
volume method were recently presented by Wang et al. cB. 
Though superficially similar to the PIC method, their 
model is approximate. They included the electronic free 
energy in the same way that we do, but no supercell is 
used for the phonon contribution. Instead they find an 
effective potential from the equation of state of the prim- 
itive, one atom unit cell, and integrate the mean field 
phonon partition function based on this effective poten- 
tial. This is a tremendous reduction in effort compared 
with the use of large supercells and finding the potential 
for displacing one atom in the supercell. They obtain im- 
pressive results for this simple model obtaining excellent 
agreement with the experimental Hugoniots, not only for 
Ta, but also Al, Cu, Mo and W. Nevertheless, it seems 
unlikely that this simply modeLwill work for lower sym- 
metry systems such as hcp-FeLD or for elastic constants. 
For example, the c/a varies with temperature in hcp-Fe, 
but the Wang et al. model would not allow for this. 

We summarize the zero pressure 300K equation of state 



parameters in Table III. The equation of state gives 
V(0,T=300 K)=124.489 bohr 3 , 2% higher than the ex- 
perimental value. Another way of looking at the discrep- 
ancy, is that the computed pressure at the room temper- 
ature experimental volume 121.8 bohr 3 is 4.1 GPa rather 



than 0. Table |IHj also shows the equation of state param- 
eters computed at the experimental volume. The main 
discrepancy is the thermal expansivity which is 35% too 
high, though this is reduced by comparing at the exper- 
imental volume. The Gruneisen parameter is similarly 
high. The origin of this discrepancy is unknown, as dis- 
cussed above, although the thermal expansivity (and thus 
7) are known to be very sensitive. Apparently this inac- 
curacy must decrease with increasing pressure, since our 
Hugoniot agrees well with experiment, up to tempera- 
tures of almost 10,000K. Perhaps our potential surface 
is not modeled accurately enough at small displacements 
and low pressures due to the very small energy differences 
involved in that regime. 



V. CONCLUSIONS 

We have studied the static and thermal equation 
of state of Ta from first principles calculations. An 
electronic topological phase transition is found around 
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200 GPa. Three different forms of thermal equation of 
state is provided as: Vinet equation of state with tem- 
perature dependent equilibrium quantities, simple linear 
temperature dependent average thermal pressure, and a 
global fit to the Helmholtz free energy F(V,T). The sim- 
ple equation of state Pth — works quite well, but 
more accuracy and insites into higher order thermoelas- 
tic parameters were obtained from the global fit in V 
and T. We find that aKr is quite constant, as has been 
seen in experiments forj-a. wide range of materials above 
the Debye temperaturescJ and has been shown for simple 
pair potentialsEa. Electronic excitations contribute sig- 
nificantly to the heat capacity temperature dependence 
of Cy , and thus to variations in the Griineisen param- 
eter 7. We find good agreement with the experimental 
Hugoniot and thermal expansivity, though the rapid in- 
crease in the thermal expansivity and heat capacity at 
high temperatures remains unexplained. 
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TABLE I. Vinet parameters for Ta isotherms. The tem- 
peratures are odd numbers since computations were done at 
multiples of k B T =0.006 Ry. 



T (K) 


Vo(T) (au) 


K (T) (GPa) 


K(T) 


E (T) (Ryd) 


0.00 


123.63 


190.9 


3.99 


-31252.3341 


947.32 


126.83 


163.6 


4.31 


-31252.3665 


2052.53 


131.38 


138.2 


4.54 


-31252.4388 


2999.85 


135.53 


120.5 


4.69 


-31252.5119 


3947.17 


140.26 


103.8 


4.83 


-31252.5925 


5052.37 


147.30 


83.1 


5.05 


-31252.6950 


5999.69 


155.44 


64.5 


5.30 


-31252.7899 


6947.01 


167.39 


45.0 


5.64 


-31252.8917 


8052.22 


192.24 


22.4 


6.25 


-31253.0210 


8999.54 


230 


9.2 


6.90 


-31253.1413 


9946.86 


289 


3.1 


7.53 


-31253.2678 



TABLE II. Global fit parameters for Ta in SI units (ex- 
cept where marked). Row is for i and column for j in Aij of 
Eq. |ll]. Note that the T = parameters are not identical to 
those in Table [l| since these were determined from a global 
fit to all results (T = and T 0) and the Table | values 
were fit to T = only. There is no practical difference for 
applications of these equations of state, within the accuracy 
of the computations. 







2.768 10 -22 -3.295 10 b 6.852 10 33 -1.751 10 63 

3.734 10" 27 -639.6 2.259 10 31 -4.524 10 58 

-1.955 10" 31 0.0247 -6.512 10 26 -1.137 10 55 

123.52 au 1.8304 10~ 29 m 3 

186.7 GPa 186.7 GPa 

4.120 4.120 

-31252.3333 Ryd -6.81261488 10~ 14 J 



1 

2 
3 

Vo 
K 
K' 
Eo 



TABLE III. Thermal equation of state parameters at ambi- 
ent conditions (300 K) (theoretical room temperature volume 
and experimental room temperature volume). 



theoretical volume experimental volume experiment 



Vo bohr 3 


124.5 


(121. 


K T0 GPa 


180 


197 




4.2 


4.07 


a lO^K" 1 


2.64 


2.38 


To 


2.09 


2.02 


Cp/R 


3.03 


3.08 




3.04ia 
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I I I I I I I I I I I I 

80 90 100 110 120 130 

Volume (au) 

FIG. 1. Static equation of state of Ta. Solid and dashed lines are LAPW GGA and LDA calculations respectively. Pseudopo- 
tential mixed-basis results are shown by the dotted line. Circles and squares are two different diamond anvil cell experiments. 
The discrepancies at high pressures may be due to strength effects in the static equation of state. 
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FIG. 2. Energy differences between the calculated and fitted data of the equation of state fits. Both Vinet and 
Birch- Murnaghan equation-of-state fits are shown. In the low pressure fits only the data with volume greater than 84 bohr 3 
were included. The extended Vinet fit included the 2 nd and 3 th pressure derivative of Ko(T). 
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FIG. 4. Vinet equation-of-state parameters as functions of temperature, (a) The minimum energy with respect to the T=0 
K value E (T = 0) = -31252.33412 Ryd, (b) the equilibrium volume, Vb(T), (c) the bulk modulus K T o{T)aX P =0, and (d) the 
pressure derivative of bulk modulus at P —Q. The vertical line shows the position of last volume point. At higher temperatures 
the parameters are Active, but still govern the high pressure equation of state. 
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FIG. 5. The thermal pressure as a function of (a) volume at different temperatures (shown in Table |l|) and (b) temperature 
from the Vinet fits. 
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FIG. 6. The thermal pressure averaged over volume as a function of temperature. Solid line is a linear fit with slope equal 
to 0.00442 GPa/K. 
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FIG. 7. Global fit to free energies, (a) Computed Helmholtz free energies (symbols) and the result of the global fit (lines) 
using parameters from Table |n|. The top curve is for K, and the bottom 9947 K. See Table [j| for exact temperatures, (b) 
Residuals of fit. Lines are labeled in order of temperature. 
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FIG. 8. The thermal expansivity as a function of (a) temperature and (b) pressure. The wide line in (a) shows the experi- 
mental zero pressure thermal expansivity (the dashed part is less well constrained) The experimental thermal expansivity shows 
an anomalous increase at high temperatures. Theory predicts the thermal expansivity to drop rapidly with pressure, and the 
temperature dependence to decrease. 
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FIG. 9. The Anderson-Griineisen parameter St as a function of (a) T, (b) P, and (c) r) = V/Vo(T = 0). The average of St 
from 0-6000 K is also shown as a thick line in (c). (d) Contours of St — K' show that this quantity is quite small for 1000-2000 
K, but increases at low and high temperatures. 
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FIG. 10. The dependence of aKr on (a) temperature and (b) pressure. The temperature dependence of cxKt changes sign 
with pressure, but the changes are small. 
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FIG. 11. Variation of the Griineisen parameter 7 with (a) temperature and (b) pressure. The temperature dependence is 
moderate. 
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FIG. 12. Variation of q as a function of (a) temperature and (b) pressure. The pressure dependence is small above 50 GPa, 
but the temperature dependence is significant at all pressures. 
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FIG. IS. Variation of the heat capacity with (a) temperature and (b) pressure. The thick curve is the experimental heat 
capacitytJ, which shows an anomalous increase at high temperatures, similar to the behavior of the thermal expansivity. 
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FIG. 14. a)Tbe line shows the computed Hugoniot for Ta (computed at the points shown as circles) The squares are the 
shock-wave datan. b) Theoretical temperatures along the Hugoniot. 
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